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Gravitational waves are generated during first-order phase transitions, either 
by turbolence or by bubble collisions. If the transition takes place at temperatures 
of the order of the electroweak scale, the frequency of these gravitational waves is 
today just within the band of the planned space interferometer LISA. We present 
a detailed analysis of the production of gravitational waves during an electroweak 
phase transition in different supersymmetric models where, contrary to the case of 
the Standard Model, the transition can be first order. We find that the stochas- 
tic background of gravitational waves generated by bubble nucleation can reach a 
S^ maximum value /io$7 gw of order (10 -10 — 10 -11 ), which is within the reach of the 

planned sensitivity of LISA, while turbolence can even produce signals at the level 
/igOgw ~ 10~ 9 . These values of /ig^gw are obtained in the regions of the param- 
eter space which can account for the generation of the baryon asymmetry at the 
electroweak scale. 
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Introduction 

Gravitational wave (GW) interferometers are presently under construction, and have as 
a possible target a stochastic background of GWs of cosmological origin (for reviews, 
see |[j, ^ [3|, U). In particular, the space interferometer LISA is expected to fly around 
2010, with an extremely good sensitivity for GWs with a frequency / between 10 -4 Hz 
and 1 Hz ||. It is quite fortunate that this is exactly the frequency range that a gravi- 
tational wave produced at a temperature T ~ 100 GeV has today, after redshifting, and 
it is therefore interesting to investigate the intensity of GW backgrounds produced at a 
possible electroweak phase transition, and to compare it with the expected sensitivity of 
LISA. 

Two main production mechanisms have been proposed for a first order phase transi- 
tion: the first is the nucleation and collision of true vacuum bubbles |6|, 0, ||, |9], [10], [U], [L2| . 



The second is the onset of turbolence as a consequence of the injection of energy into the 
primordial plasma [|T^, [T3|] . 

These mechanisms are only effective for a strongly first order phase transition. For 
the Standard Model, the strength of the electroweak transition depends of the Higgs 
and W masses, and can be investigate perturbatively only if itlh <C mw- The physical 
situation m# > mw requires instead non-perturbative lattice computations, which have 



shown |L4| that unfortunately there is no phase transition at all, but rather a smooth 
crossover. Therefore, in the Standard Model, no GWs are produced by this mechanism. 

The strength of phase transitions, however, is strongly model- and parameter-dependent 
(for review, see JTJI), and it is therefore meaningful to investigate whether in extensions 
of the Standard Model one can have a background of interest for LISA. In this paper we 
investigate the production of GWs in supersymmetric extensions of the SM. The strength 
of the phase transition has been investigated in connection with the generation of the 
baryon asymmetry [|l(|, and it is known that in the Minimal Supersymmetric Standard 



Model (MSSM) a strong enough phase transition requires light Higgs and stop eigen- 
states [[17]]. If the Higgs is heavier than about 115 GeV, stronger constraints are imposed 
on the space of supersymmetric parameters. If one goes beyond the MSSM, introducing 
an additional gauge singlet in the Higgs sector (the so-called next-to-minimal SSM), it 
is known that the strength of the transition can be further enhanced ||18|| . 

The paper is organized as follows. In section [1] we review the physical processes that 
can lead to a stochastic background of gravitational waves during a phase transition; in 
particular, in section [L3| we describe the strategy of our computation. Sections II and 
III are dedicated to the analysis of GW production, from turbolence and from bubble 
collisions, in two specific supersymmetric models: the Minimal Supersymmetric Standard 
Model (MSSM) and the Next-to-Minimal Supersymmetric Standard Model (NMSSM); we 



present a detailed analysis that extends results already reported in |]19| . Our conclusions 
are presented in section IV. 



1 GW production at phase transition 

1.1 GW background from bubble collision 

In first order phase transitions the Universe finds itself in a metastable state, which is 
separated from the true vacuum by a barrier in the potential of the order parameter, 
usually a scalar field 0. True vacuum bubbles are then nucleated via quantum tunneling. 
In strongly first order phase transition the subsequent bubble dynamics is relatively sim- 
ple: once the bubbles are nucleated, if they are smaller than a critical size their volume 
energy cannot overcome the shrinking effect of the surface tension, and they disappear. 
However, as the temperature drops below a critical temperature T c , it becomes possi- 
ble to nucleate bubbles that are larger than this critical size. These 'critical bubbles' 
start expanding, until their wall move at a speed close to the speed of light. The energy 
gained in the transition from the metastable state to the ground state is transferred to 
the kinetic energy of the bubble wall. As the bubble expands, more and more regions of 
space convert to the ground state, and the wall becomes more and more energetic. At 
the same time, it also becomes thinner, and therefore the energy density stored in the 
wall increases very fast. As long as we have a single spherical bubble, this large energy of 
course cannot be converted in GWs. But when two bubbles collide spherical symmetry 
is broken and we have the condition for the liberation of a large amount of energy into 
GWs. In particular, there are two possible 'combustion' modes for the two bubbles: det- 
onation P0| , that basically takes place when the boundaries propagate faster than the 
speed of sound, and deflagration, when instead they move slower. In the first case there 
is a large production of GWs 0, ||, |12| . The dynamics of a particular phase transition 
is governed by the effective potential for the scalar field(s) driving the transition itself. 
Given this effective potential, two are the basic quantities which play a role in the de- 
termination of the GW background generated during a first-order phase transition, and 
in which all the relevant features of the potential are encoded. The parameter a gives a 
measure of the jump in the energy density experienced by the order parameter </> during 
the transition from the false to the true vacuum and it is the ratio between the false vac- 
uum energy density and the energy density of the radiation at the transition temperature 
T*. The parameter f3 characterizes the bubble nucleation rate per unit volume, which can 
be expressed as Y = r exp(/3t) flO|| . Thus /3" 1 = T/T represents roughly the duration 
of the phase transition and the characteristic frequency of the gravitational radiation at 
time of production is expected to be 2irf ~ (3. The stronger is the transition, the larger 
is a and the smaller is (3 because of the larger amount of supercooling experienced by 
the system before the transition. 

Once one has determined a and /3, the gravitational radiation resulting from bubbles 
collisions can be calculated, for a strongly first order transition and for detonation com- 
bustion mode, in function of these two parameter only. We refer the reader to Ref. [T2 



for a detailed derivation and we summarize here only the strategy. The stress-energy ten- 



sor for a bubble expanding as a detonation front at velocity v and of size R ~ (3~ l v takes 
a simple form, its relevant part coming from the wall's contribution only. Its amplitude 
is proportional only to the kinetic energy density of the wall, which is a fraction k of the 
energy gained in the transition, namely the false vacuum energy density e times the bub- 
ble's volume R? ~ (/3~ l v) 3 . Then the energy radiated during the collision of two bubbles 
is calculated from their stress-energy tensors using the envelope approximations, while 
the theory of relativistic combustion supported by numerical calculations gives k and v 
as a function of a. Taking in account the redshift from time of production up to now, 
and evaluating numerically the contribution from many randomly nucleated bubbles in 
a sample volume, one finally obtains the energy density of the radiation produced. It is 
convenient to express it in terms of the dimensionless quantity 

pc a log/ 

where p^ is the energy density associated to GWs, / is their frequency and p c is the 
present value of the critical energy density, p c = 3Hq / (8ttG n) , with H = 100h Km/sec 
Mpc; h parametrizes the uncertainty in H (note that the combination h,QQ gw is inde- 
pendent of ho). LISA is expected to reach a sensitivity of order 

/igftgw ~ 1(T 12 at / = 1 mHz . (2) 

At this frequency a cosmological signal could be masked by an astrophysical background 
due to unresolved compact white dwarf binaries. Its strength is uncertain, since it de- 
pends on the rate of white dwarf mergers and it is estimated to be || 

/igtyjw ~ kt 11 . (3) 

At a frequency / ~ 10 mHz the LISA sensitivity is expected to be of order /iQ^ gw ~ 10 -11 , 
and the astrophysical background is expected to be below this value. 

At the peak frequency /peak, the intensity h^Q^ of the radiation produced in bubble 
collisions is given, in terms of the parameters a and (3 discussed above, by |L2] 

sw U 1, l + 0.7a 1 + a) \0) \0.2B + v*)\gJ ' [) 

where H* and g* are respectively the Hubble parameter and the number of relativistic 
degrees of freedom at the time of transition and v is the velocity at which the bubble 
expands. In the following we will use the value of the velocity v = v(a) as given in Ref. 
[ [T2]| for bubble detonation. 

The detailed computations of a and (3 will be the main issue of this paper. However, 
to have a first idea of the numbers involved, we note that early analytical estimates by 
Hogan [^l], indicated that it is difficult to obtain values of H*/j3 larger that O(10~ 2 ) 



(this estimate will be confirmed by our numerical results) and therefore, even for large 



values of a, one should expect a signal at most of order /igfi gw ~ 10~ 10 . This maximum 
value can be obtained only for strongly first order phase transition. The weaker the 
transition, the smaller are both H*/f3 and a, and therefore we have further suppressions. 
It is clear that, since the value /ij^gw ~ 1CT 10 is already close to the limiting sensitivity 
of LISA, one can hope to find an observable signal only in some specific models, and 
in special regions of the parameter space; an observable signal will not be a generic 
property of extensions of the Standard Model. However, we know that the condition 
for strongly first order phase transition is the same that is required for producing the 
baryon asymmetry at the weak scale. Therefore the models and the regions of parameter 
space which produce an interesting GW signal, albeit quite specific, are also the most 
interesting if the baryon asymmetry has been indeed generated at the weak scale, and 
it makes sense to perform a detailed scanning of the parameter space in search of the 
strongest GW signal. 

The present (i.e. properly red-shifted) peak frequency of the GW background is 



determined by f3, by .£/*, g*, and by the transition temperature T* [O 



JL\ ( T - \ (Ji.Y" 



At frequencies lower than / pea k the energy density increases like f 25 , while above drops 
off more slowly |TT| . As it will be shown in the following, typical values for /3/H* for the 



electroweak phase transition are between 10 2 and a few times 10 3 , with T* ~ 100 GeV. 
This gives a frequency / pea k i n the range (10~ 4 — 5 x 10~ 3 ) Hz, which is precisely the 
range in which LISA achieves its maximum sensitivity. 

1.2 GW background from turbolence 

During a first order transition there is another very general and possibly powerful source 
of gravitational waves. Part of the energy gained in the transition from the metastable 
state to the ground state is used to heat up the plasma, and another fraction is converted 
into bulk motion of the fluid. If the Reynolds number of the Universe at the phase 
transition is large enough, then this results in the onset of turbolence in the plasma, 
and consequently in the production of GWs. To compute accurately the amount of 
gravitational waves pruduced by turbolence is certainly a very difficult task. In the 
following, we will use the estimates for the characteristic frequency and for the intensity 



given in ref. (12 



/u - 2.0 x 10-" Hz — ( — I ( T * ) f-^V /l ' (6) 

J .. I □■ I V 100 GeV/ VIOO/ V ; 




^^lO^i^Jvvt^y* . (7) 

Here v® is the fluid velocity on the largest length scales on which the turbolence is being 
driven, and v is again the velocity of the bubble wall. For a detonation, w can be 



estimated to be vq ~ ^[kol for a weak transition (small a) and v o ~ 1 for a strong 
transition fL2f . Strictly speaking, eq. (^) is derived in the regime of nonrelativistic fluid 
velocities (small v o) • A more detailed analysis of the GW production by turbolence has 
been recently announced as in preparation |13[ . 

Substituting again typical (3/H* values of order 10 2 — 10 3 , and assuming large enough 
Vq (detonation), we obtain a peak frequency around 10 -4 — 5 x 10~ 4 Hz, again well within 
the window of sensitivity of LISA. This is therefore another, very interesting, mechanism. 



1.3 Computational method 

Before analysing specific models, we describe how one can compute the relevant quantities 
appearing in Eqs. @) and ©, namely T*, a and f3/H*. We will perform a perturbative 
computation based on the thermal effective potential. This method is known to receive 
large non-perturbative corrections when the phase transition is radiatively induced |22|, 
23| , p4| |25| , |26|| . In this case our results are really an upper bound on the GW production. 
However, the only situation in which we will find a signal compatible with LISA is when 
the potential has a metastable minimum already at tree level, and in this case our results 
are not expected to get important non-perturbative modifications. 

The typical temperature dependence of a potential V((p,T) for the Higgs scalar 
field(s), for a first-order transition, is shown in Fig. [I] (the figure refers for simplicity 
to one-dimensional potential, but the following discussion is valid also if several scalar 
fields are involved). 
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Figure 1: Typical temperature-dependence of the potential V for a scalar field driving 
a first order transition. The dotted curve refers to very high temperature, while longer 
dashed lines refer to lower temperatures. T\ is the temperature at which a minimum at 
(4>) 7^ develops; at T dcg the two minima are degenerate; at T dcst the origin ((f)) = 
becomes unstable. The actual transition temperature T* is between 7d est and T^ e „. 
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In the case of a first order transition the rate of tunneling per unit volume T from the 
metastable minimum to the stable one is in general suppressed by the exponential of an 
effective action, T = T e~ s with 



S 



drd x 



l(f+iw)' + w 



(8) 



where r is the euclidean time. In the case of field theory at finite temperture T, one is led 
to consider Euclidean field theory periodic in imaginary time with period T _1 ; moreover, 
in a cosmological context the temperature T of the Universe decreases in time. Thus in 
our case the action (H) acquires a dependence on cosmic time t, and has to be computed 
in the space of functions periodic in Euclidean time r, 



S(t)= / drd 3 x 
Jo 



Hi +w+mT) 



(9) 



where V(4>, T) is the effective potential, shifted in such a way that V(0, T) = 0. One 
then expands S(t) around the transition time £*, S(t) ~ S(t*) —/?(£ — £*) and identifies (5 
with —{dS/dt)t„- From V ex exp(— S(t)) it then follows that T ex exp(/3t). In a radiation 
dominated Universe this gives 



TT 



d(S 3 /T) 



dT 



(10) 



r i\ 



where we have considered the large T limit of 
Euclidean action 



and have defined S3 as the spatial 



Ss(T) 



d 3 x 



|(V0 6 ) 2 + U(0 fe ,T) 



(11) 



computed for the configuration of the scalar field(s) </>& describing the bubble. The rate 
per unit volume of nucleation of a critical bubble is therefore 



T = r exp 



T 



(12) 



where To is of the order of T 4 for a finite temperature field theory. 
The Euclidean action for a spherical configuration is 



S 3 (T) = 4tt f drr 2 



d( 



dr 



+ V(<p b ,T) 



the bubble configuration is thus the solution of the equation(s) 

d 2 6 h 2 dd>h dV 



h 2d4b 

dr 2 r dr 







(13) 
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Figure 2: Typical profile of a bubble solution 0&(r) interpolating from the true vacuum 
phase at r = to the false vacuum at r = oo. The solution shown corresponds to 
the critical bubble at the transition temperature in the MSSM for sin 2 /3 M gsM = 0.8, 
^Higgs = HO GeV and m st0 p = 140 GeV, see section II. 



supplemented by the boundary conditions 



dr 











(15) 



namely true vacuum inside the bubble and metastable vacuum outside. A typical bubble 
profile solution is shown in Fig. |[ In order to find the "escape point" (i.e. the value 0(r = 
0) = e , implicitily determined by the two conditions (|i~5"D) one uses the "overshooting- 
undershooting" method illustrated in fig. §. Equation (|14]) is the classical equation of 
motion for a point particle subject to a potential — V(<f>) and to a velocity- and time- 
dependent friction force, if 0(r) stands for the trajectory of the particle. If the particle 
starts at rest precisely from the escape point, it will have just enough energy to overcome 
the friction force and to come at rest at <p — at infinite r (case (a)). If instead it starts 
at the right of e (overshooting), it will continue toward = — oo (case (6)). If it starts 
at the left of </> e (case (c), undershooting) it will experience damped oscillation around 
the minimum of the inverted potential. Thus one determines the escape point by trials 
and errors, lowering the starting value if one gets a solution of type (6), and increasing it 
if one gets a solution of type (c). Once the escape point is found, the extremal action ([13]) 
can be computed for the corresponding solution. The bubble solution only exists in the 
range T dest < T < T deg , see fig. []]; for a one- dimensional potential the quantity S 3 (T)/T 
behaves as shown in fig. f|: at T > Td eg the transition cannot take place, because it is not 
energetically favorable to go from one minimum to the other, thus £3 — >• oo and r — > 0. 
Conversely, for T approaching Td es t the transition becomes more and more convenient, 
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Figure 3: Overshooting-undershooting method. 
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Figure 4: The plot shows the rapid variation experienced by the action S 3 (T)/T (evalu- 
ated on critical bubble solutions) in a small range of temperature: S < 3(T)/T moves from 
to oo going from T = T dcst to T = T deg , with (T deg - T dest )/T dest ~ 1%. The par- 
ticular case shown in the figure refers to the electroweak transtion in the MSSM with 
"^Higgs = HO GeV, m st0 p = 145 GeV, sin 2 /3MssM = 0.8, see section II. 

S3 — > and r — » r . The temperature T* at which the transition takes place is computed 
by comparing the probability of bubble nucleation per unit time and unit volume ([12]) 
with the expansion rate of the universe at the corresponding temperature. The condition 
that must be satisfied is that the probability for a single bubble to be nucleated within 
one horizon volume is 0(1), 

e -S 3 (T)/T „ 1 (16) 




For the electroweak transition, in which characteristic temperatures are O(100 GeV), the 
condition above is well approximated by S3 (T*)/T* ~ 140. 

Once one has found T*, the parameter a is readily computed from the definition given 
above. Calling v(T) the vacuum expectatuion value of the Higgs field at the true vacuum 
at temperature T, then the quantity 



e*= -V(v(T),T) + T-^V(v(T) } T) 



(17) 



is the vacuum energy (latent heat) density associated with the transition (remember that 
the effective potential V((f>,T) is the free energy density). Thus a is the ratio between 
this and the energy density of the radiation, 

a = ^vr (18) 

10 



while (3/H* is given by eq. (|T0|) . Notice that in order to compute a it is sufficient to 
know the transition temperature T* and the potential V(<p,T), while the determination 
of /3/H* is possible only if one knows how the action Ss(T) computed on the exact bubble 
solutions depends on T in a neighbourhood of T*. 

2 The Minimal SuperSymmetric Standard Model 

As we have already noticed, extensions of the SM are required to obtain a first-order phase 
transition at the electroweak scale, because in the SM, in the physical case m H igg S > mw, 
there is no phase transition. 

The Higgs sector of the MSSM requires two complex Higgs doublets, with opposite 
hypercharges 

' #M u _ ( h.; 



*-UJv *-UV (19) 



The tree-level potential is 



V® = mJHl^+mlHl^+mK^^+h.c.^g 2 (H ] 2 aH 2 + HlaH^+^g' 2 (h\H 2 - H\H x 

(20) 
where the products between doublets are the usual SU(2) invariant products, the <Tj are 
the Pauli matrices, g and g' are the gauge couplings and the m 2 are real parameters, not 
necessarily positive. In fact, in order to have the gauge group SU(2)xU(l) broken at 
zero temperature, one needs the quadratic term to be negative along some directions in 
the Higgs fields space. 

We define V\ = (H®), vi = (H 2 ) and tan/^MssM = vijv\. In the following we will 
indicate the angle /9mssm simply by (3: it should not be confused with the parameter /3 
discussed in the previous sections that sets the time scale of phase transition. 
As usual, the W mass fixes v 2 = v\ + v\ to the value (246 GeV) 2 . If the mass m^ 
of the CP-odd field in the Higgs sector is much larger than mw, only one light Higgs 
scalar is left and its tree-level potential is identical to the one in the SM, namely 
Vq{4>) = — ^-0 2 + |0 4 . This will be the case considered in this section. When m^ ~ m w 
the full two-Higgs potential should be considered, but the strength of the phase transition 
is weakened ]l7j . As we will see in the following, the one-loop thermal corrections to the 
effective potential make the quadratic term positive at high temperature and create 
a negative cubic term due to loops of the massive bosons in the theory. Due to the 
presence of this cubic term and the positivity of the quadratic one, there exists a range 
of temperature in which the point <ft = is a local minimum separated from the true 
simmetry-breaking one by a small potential barrier, that is precisely the set-up for a first- 
order phase transition. The strength of the phase transition is enhanced by the presence 
of new bosons coupled to the Higgs, a significant role being played by the right-handed 
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stop £r, which is - apart from the Higgs itself - the lightest scalar in the theory and has 
the largest Yukawa coupling to the Higgs 0. 

To study the amount of gravitational waves generated during the electroweak phase 
transition within the MSSM we have made use of the thermal potential corrected up to 



two- loop level |29|]. However, as discussed above, in this radiatively induced case, non- 
perturbative effects are also important. The results of ref. |2|| indicate that the two-loop 
perturbative computation overestimates the amount of supercooling, so that the results 
of this section are really upper bounds on the GW production. 

Following |2£| , we summarize here the one and two- loop finite temperature contribu- 
tions to the effective potential for the Higgs field 0. At one loop, the dominant corrections 
come from the gauge bosons Z and W, the top quark and its supersymmetric partners. 
We will work in the limit in which the left handed stop is heavy, ttiq £ 500 GeV; lower 
values of rriQ make the phase transition stronger, but give large contribution to the p- 
parameter. For right-handed stop masses below (or of order of) the top quark mass, and 
for large values of the CP-odd Higgs mass, m^ ^> rn Z) the one- loop effective potential 
admits the high temperature expansion [29| 

3/2" 



W) + W,T) 



m 2 (T) 



2 A(T) 
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-T 



E, 



SM 



+ (2iv c ; 



(mL p (0)+n stop (T) 



12vr 



(21) 

where the various quantities are defined as follows: N c = 3 is the number of colors; E$m 
is the cubic term coefficient in the Standard Model case, 



E, 



SM 



1 

3 



'2771^ + m|' 
2wv 3 



(22) 



with my/ and mz standing for the physical masses of the gauge bosons at zero temper- 
ature. Here m stop (0) is the lightest stop mass, which depends explicitely on the VEV of 
the Higgs (ft and is approximately given by 



m, 



stop 



(<p) ^mfj + 



m\ 



nr 



0.15-f cos 2/3 + — f ( 1 \ 



A\ 



m. 



(23) 



where m t is the zero-temperature top quark mass, mfj is a parameter of the model - the 
soft mass squared, not necessarily positive - and A t is the stop mixing parameter (which 
for simplicity will be set to zero in our numerical analysis). Here n stop (T) is the finite 
temperature contribution to the right-handed stop self-energy, 



n st op(T) 



1^ + \hi 



1 + sin' 



(3 (l - Al/ml) 



T 2 + -- 



cos2(3\) g"T\ (24) 



,3 U 

where g s is the strong gauge coupling and ht is the top Yukawa coupling; in our notations 
m t = h t vsmj3/\/2. Finally m 2 {T) is given by 



2 m Higgs + 



T 2 



-m 



Higgs 



+ -ttvL + —m\ r + m 



6 



i'W 



12 
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(25) 
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and A(T) ~ ' m Higgs/' i;2 apart from logarithmic corrections. As was observed in Ref. |30 



the phase transition strength is maximized for values of the soft breaking parameter 



*l 



-il stop (T), for which the coefficient of the cubic term in the effective potential, 



E 



/if sin 3 /? f 1 - A 2 t /m 2 Y /2 



which governs the strength of the phase transition (in this case v(T*)/T* ~ 4E/X), can 
be one order of magnitude larger than E$m POR - However, it was also noticed that 
such large negative values of mfj may induce the presence of color breaking minima 
because the effective stop mass at finite temperature m 2 top ((p) + n stop (T) may become 
negative Q, |3~0|. The problem arises for small values of the VEV of </> and of the 



temperature T, as one can easily infer from Eqs. (^) and (|24|). Thus, the check one 
has to perform is that for T > T* - i.e. before the phase transition takes place - the 
sum m 2 top (((f))) + n stop (T) = m 2 ; + n stop (T) is positive, the equality being motivated 
by the fact that (</>) = for T > T*. Once the above condition is satisfied, the color 
breaking minima problem is absent for T < T* too, because the VEV (<f>) gets a large 
value 0(100 GeV) after the phase transition. 

Before performing the full computation with the two-loop corrections to the potential, 



we proceed to a rough computation of the quantity a defined in section |L1| , using the 
one-loop effective potential fl2T|), in order to understand the dependence of the strength of 
the phase transition on the various parameters of the model, and to identify the regions 
of the parameter space where it is appropriate to perform more detailed computations. 
Under the assumptions mentioned above, the most relevant parameters in the game 
are the Higgs mass mHi ggs , the right-handed stop mass m stop and the zero temperature 
ratio between the vacuum expectation values of the two neutral Higgses tan /?. However, 
it is well known that the strength of the transition has a very slight dependence on 
tan/?, with slightly stronger transition at high (i.e. greater than about 0.6) values of 
sin 2 /?. In our computations we will set both sin 2 /? = 0.8 and sin 2 /? = 0.3, obtaining 
very similar results. Our strategy has been the following. For every value of mniggs and 
"^stop we computed the temperature Td es t at which the origin = gets destabilized; 



then we computed a setting T* = T dest in the definition ([TJ). The actual temperature 
T* is in the range between the "destabilization temperature" Td es t and the "degeneracy 
temperature" Td eg , as discussed in section [T73| . For the model under study this range is 
very small, (T deg — T dest )/T dcst ~ (1 — 3) %. So even taking T* to be equal to T dest leads 
nevertheless to the right order of magnitude in the calculation of the parameter a. Our 
results are summarized in figs. [| and |], which show contour plots for a as a function of 
mHiggs and m stop , together with the region forbidden by the condition of absence of color 
breaking minima discussed above. From figs. |5] and |6| we see that in order to have large 
values for a one needs small Higgs and stop masses; this agrees with the well known fact 
that the strength of the electroweak transition is enhanced by taking light Higgs and 
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Figure 5: Contour plot of a as a function of the Higgs and stop masses for sin 2 (3 = 0.8. 
The shaded region is forbidden by the request of absence of color breaking minima, a is 
computed at the destabilization temperature, as discussed in the text. 
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Figure 6: Contour plot of a as a function of the Higgs and stop masses for sin 2 (3 = 0.3 . 



stop. We see also that the strength of the transition is very slightly dependent on sin 2 /3; 
thus for defmiteness we fix sin 2 (3 = 0.8. 

The most important two loop corrections are of the form cf 2 log(0) and are induced 
by the Standard Model weak gauge bosons as well as by the stop and gluon loops |ER [TT| 

2 / 



V 2 (<f>,T) 



L2rp2 



T 2 



32vr 2 



^ 2 -3/4sin 4 /3(l-4 
16 \ rrin 



+ 8g 2 s h 2 t sin 2 1- 



log 



A 



H 



(27) 

where the first term comes from the Standard Model gauge boson-loop contributions, 
while the second and third terms come from the light supersymmetric particle loop 



contributions. The scale A# depends on the finite corrections p8[ and is of order of 
100 GeV: given the slight logarithmic dependence of V 2 on A#, in the following we will 
set A H = 100 GeV. The complete potential we use is thus V((f), T) = V ((p) + Vi((f>, T) + 

Next, having fig. [^ in mind, we proceeded to an accurate numerical computation of 
the parameters a and (3/H* characterizing the strength of the phase transition: for any 
given choice of the masses mHi ggs and m stop , we have first numerically computed the 
nucleation temperature T* by imposing that, for the Higgs field configuration describing 
the nucleated bubble, the condition Sz{T*)/T* ~ 140 is satisfied. Then we have computed 
the parameters a and (3/H* through Eqs. dl8|) e (0). Our results are summarized in 
Figs. |7| and ||. 
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Figure 7: Plot for a (dashed line), (3/H* (dotted line), and /ijj^gw from bubble collisions 
(solid line), as functions of the stop mass, for a 110 GeV Higgs and sin 2 /3 = 0.8. 

The intensity of the produced gravity waves rapidly grows as the masses decrease. 
This is explained by the fact that HqQ sw scales roughly as a 3 //3 2 and a linear increase of 
a and decrease of f3 leads to a rapid power-law growth of the intensity. Unfortunately, 
the general prediction is that the intensity of the gravitational waves produced during 
the MSSM phase transition is too small for LISA. For instance, taking a Higgs mass of 
110 GeV, the right-handed stop mass of 140 GeV and sin 2 /? = 0.8 (see fig. |?D, we find 
a ~ 3 x 10 -2 and (3/H*~ 4 x 10 3 , leading to h^Vt^ ~2x 10~ 18 . Notice that one is not 
allowed to lower too much the right-handed stop mass because of the problem of color 
breaking minima described above. 
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Figure 8: Plot for a (dashed line), (3/H* (dotted line), and h^Q^ from bubble collisions 
(solid line), as functions of the Higgs mass, for a 155 GeV stop and sin 2 /3 = 0.8. 

We have also estimated what happens if we lower the Higgs mass down to (the already 
excluded value of) 80 GeV, setting the right-handed stop mass at 155 GeV - which is 
the lower value compatible with the absence of color breaking minima - and sin 2 (3 = 0.8, 
see fig. ||. We obtain a ~ 0.1 and /3/i7*~ 2 x 10 3 , giving h^Q^ ~2x 10~ 16 , a signal 
still not relevant. The situation does not improve when both Higgses are involved in 
the transition because the strength of the phase transition is weaker. An uncertainty in 
this estimate is due to the determination of v, the velocity of expansion of the bubble 



discussed in section 1.1. If the phase transition is not strong enough, then v is subsonic, 



so that the value of hf,VL PVJ is further suppressed. 



tJ^gw 
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It is now easy to estimate also the amount of GWs produced by turbolence at the 
electroweak transition, in the MSSM case. Substituting in eq. (|^) tha values of a and 
P/H* found above, we obtain the plots shown in figs. ^| and [U]. Again, the results are 
too small to be of interest for LISA. 
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Figure 9: hlQ gw from turbolence, as a function of the stop mass for a 110 GeV Higgs 
and sin 2 /3 = 0.8. 
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Figure 10: /ig^gw from turbolence, as a function of the Higgs mass for a 155 GeV stop 
and sin 2 (3 = 0.8. 



3 The Next-to-Minimal Supersymmetric Standard 
Model 
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3.1 The model 

The situation improves considerably if we enlarge the MSSM sector adding a complex 
gauge singlet iV |?7|. This is the so-called Next-to-Minimal Supersymmetric Standard 
Model (NMSSM) and is a particularly attractive model to explain the observed baryon 
asymmetry at the electroweak phase transition. The relevant part of the superpotential 
is given by 

W = XH t H 2 N - -N 3 , 

where now the supersymmetric //-parameter of the MSSM is substituted by the the 
combination \(N), and k is a free parameter. The corresponding tree level Higgs potential 
reads V = Vf + Vd + V 80 ft, where 



V F = \X 



2 



WdHtf + \H 2 \ 2 ) + \H X H 2 \ 2 + k 2 \N\* - (Xk*H x H 2 N 2 * + h.c.) (28) 



Vd = 9 -^(\H 2 \ 2 -\H 1 \ 2 ) 2 + 9 -^\hIh 2 \ 2 (29) 

o Z 

V soit = ml 1 \H 1 \ 2 +ml 2 \H 2 \ 2 + m%\N\ 2 -(\A x H 1 H 2 N-^kA k N 3 + h.c^ .{30) 

The presence of the cubic supersymmetry breaking soft terms proportional to the param- 
eters A\ and A), already at zero temperature makes it clear that within the NMSSM it is 
quite easy to get a very strong first-order phase transition at the electroweak scale |TJ[ . 



The order of the transition is determined by these trilinear soft terms rather than by the 
cubic term appearing in the finite temperature one-loop corrections and the preservation 
of baryon asymmetry after the phase transition is possible for masses of the lightest scalar 
up to about 170 GeV. At the same time, since the transition is induced by the tree-level 
terms rather than by radiative corrections, non-perturbative effects are non expected to 
change dramatically the results of the perturbative computation. 

By a redefinition of the global phases of Hi and N it is always possible to take 
\A\ and kA k real and positive, while by an SU(2)xU(l) global rotation one can put 
v^ = (H^) = and v 2 = (H®) G R + . Assuming CP conservation, also v® = (Hi) and 
x = (N) are real. 

Before introducing finite temperature corrections, we summarize here some aspects 
of the tree level potential. First, it contains seven free parameters, but imposing the 
stationarity conditions in (H^H^N) = (vi,v 2 ,x) with the constraint vf + v\ = v 2 = 
(246 GeV) 2 , we can express the soft masses in terms of the six parameters A, k, A\, Ak, 
tan/? and x, 



J 2 + a 12 
m 2 Hl = X(A X + kx)x tan (3 - X 2 (x 2 + v 2 sin 2 (3) - - / v 2 cos 2(3 



4 

,2 i „/2 



m 2 H2 = \(A X + kx) x cot&n(3 - X 2 (x 2 + v 2 cos 2 (3) + - A 9 v 2 cos 2(3 (31) 



4 



v 2 

,2 _ \ A U •„ OP , U A ™ \2„.2 o7,2„2 , \7_.2, 



m z = \A X — sin 2(3 + kA k x - AV - 2k V + Xkv 2 sin 2(3 
2x 
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The above equations do not guarantee that (v±,V2,x) is the global minimum of the tree 
level potential: for each choice of the parameters we must verify that this is indeed the 
case. 

Concerning the physical masses of the 6 physical Higgs fields, one has to consider 
the matrix of the second derivatives of the tree-level potential. In order to avoid the 
formation of dangerous directions of instability in the space of charged and pseudoscalar 
Higgses, in the analysis of the strength of the transition for each choice of the parameters 
one has to verify the positiveness of their mass matrixes. 

Finally, barring the possibility that the transition occurs along CP-violating direc- 
tions, the potential becomes a function of three real scalar fields only, (KeH °, ReH®, ReiV) 
which in the following will be denoted by (</>i, </> 2 , ^3). 

In our numerical analysis we made use of the tree- level potential Vp + Vd + V so f t , 
restricted to the three neutral scalar fields left, plus the one- loop corrections appearing 
at finite temperature. Overall, we have six free parameters: the coupling parameters A 
and k; the soft-breaking mass terms A\ and A^, the zero-temperature vacuum expectation 
value of the singlet x and tan/?. The thermal corrections to the potential involve a linear 
term, 



and a bilinear one, 



Vsin(2/3)AA A 3 , (32) 

o 



^ (C Hl <Pi + C H2 4>l + C N <% + 2C 12 fafa) , (33) 



where the C"s are constants depending on A, k and tan/? [IS 



3.2 The landscape of minima 

Since we have now to deal with three scalars and six free parameters, the NMSSM case 
is considerably more complicated than the case of one light Higgs in the MSSM. Besides 
increased computational difficulty, new features arise. 

First of all, the strength of the transition is dominated by the tree-level trilinear term 
rather than by loop corrections, and therefore it is not anymore directly related to the 
Higgses and stop masses (neither is related in a simple way to the six parameters of the 
tree level potential). 

More in general, the landscape of local and global minima in this large parameter 
space is more complicated. Indeed, for some choices of the parameters the destabiliza- 
tion does not even take place, while for other values of the parameters, at T = Td es t 
the origin gets destabilized along a direction which does not connect the origin to the 
true vacuum, but to a local minimum, separated from the global one by a barrier: in 
these cases the system first rolls down to the new "shifted" false vacuum (second order 
transition), and next it tunnels through the barrier which separates it from the true vac- 
uum (first order transition). These two minima are separated by a potential barrier in 
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the three-dimensional field space (0i, 02,03); tunneling throughout the barrier will take 
place through the trajectory in this field space which leads to the least Euclidean action 
( |13D and this in general will not be a straight line. Our approximation from now on 
is that this least-action trajectory is in fact a straight line joining the two minima: in 
other words, after having exactly computed the two minima, we reduce the problem to 
a one dimensional system along the direction joining these two points, thus simplifying 
the computation of bubble solutions. This perhaps leads to a small overestimate of the 
strength of the transition. 

In our analysis we have focussed on those regions of the parameter space which pre- 
vious studies have shown to give rise to a large baryon asymmetry; our strategy has been 
the following. First, we made a quicker scanning of wide regions in the six dimensional 
parameters space. At this stage, it would be highly impractical to reconstruct the bubble 
profile and compute a and (3 for each value of the parameters explored; instead, in order 
to have an idea of what are the most interesting regions, we computed, at the degener- 
acy temperature and (if it exists) at the destabilization temperature, the height of the 
potential barrier and the energy density difference between the false vacuum (possibly 
shifted, as discussed above) and the true vacuum. 

These first estimates allowed us to select the regions in the parameters space which 
present a strong phase transition; in particular, it turns out that the most important 
parameters in the game are the soft breaking mass terms A\ and Aj.- We then moved to 
a more accurate analysis restricted to some of the interesting regions found. 

The situation is now much more complicated than in the MSSM case, becuase of 
the higher dimensionality of the Higgs space. Furthermore, now the barrier is much 
higher than in the MSSM, and therefore the system has to experience an higher degree 
of supercooling before the transition can take place. Thus, while in the MSSM the 
separations between the temperature of degeneracy for the two vacua Td eg , the actual 
transition temperature T*, and the origin-destabilization temperature Td es t, are of order of 
a few GeVs, in the NMSSM these intervals can be sensibly larger (and indeed, as discussed 
above, for some ranges of the parameters the destabilization can even be absent). This 
fact makes it harder the compute precisely the value transition temperature T*, and thus 
of a. In order to get a contour plot for the function a(A\, A k ) we proceeded as follows. 
The transition takes place when S^(T)/T ~ 140, S3 being the sum of a kinetic term and 
a potential one (see (0)). These two terms are likely to be of the same order if evaluated 
on the bubble-like solution of the equations of motion, and thus Ss(T) oc Jdrr 2 V(4>b, T). 
The integral of the potential extended to the bubble receives a contribution mainly from 
the wall proportional to the maximum V max of the potential and a contribution from 
the inside of the bubble proportional to the minimum V m { n , i.e. Ss(T) ~ A ■ V max (T) + 
B ■ V min (T), with A and B depending on the exact solution b (r). The condition to 
determine T* thus becomes A ■ V max (T*) + B ■ V min (T^) ~ 140 T*. We assumed constant 
A and B in each small region we explored, thus fitting their value with several exact 
solutions 0b(r). As in the MSSM case, once this (approximate) transition temperature 
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Figure 11: Rough contour plot for a as a function of A x and A k , for x 
A = 0.83, k = 0.67 and tan/5 = 2 . 



520 



350 GeV, 



T* is evaluated, one easily computes a from its definition (|TB|). A result of this procedure 
is summarized in fig. [XT], which shows a rapid growth of a versus A\ and A/., with x = 350 
GeV, A = 0.83, k = 0.67 and tan/3 = 2. Next we moved to computing exactly T*, a 



and (3/H* for specific values of A\ and Aj~ in the region of fig. |TT|, again by numerically 
computing the temperature for which the bubble solution 4>b{f) gives S < 3(T)/T ~ 140 
and then applying relations (|T8|) and (0). Typical results are shown in Figs. [L2| and[TJ], 



fc- 



where we plot a, (3/H* and /igfi gw as functions of A\ and A 

Note that, since the transition in the NMSSM is strongly first order, it is correct to 
use the expression for v (a) computed for detonation bubbles [Q. The results shown in 
the figures are not peculiar of the values of the parameters chosen: similar results can be 



found in completely different regions of the parameters space, see for instance fig. [16 



An example of a region of the parameter space in which the transition happens 
through the nucleation of bubbles subsequent to a smooth roll-down is shown in fig. Il8l 
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Figure 12: Plot for a (dashed line), (3/H* (dotted line), and hlfl gw (solid line) from 
bubble collisions, as functions of A\ for Ak = 480 GeV, x = 350 GeV, A = 0.83, k = 0.67 
and tan/3 = 2. 
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Figure 13: tiffin from turbolence as a function of A\ for A k = 480 GeV, x = 350 GeV, 
A = 0.83, k = 0.67 and tan/3 = 2. 
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Figure 14: Plot for a (dashed line), (3/H* (dotted line), and hlfl gw (solid line) from 
bubble collisions, as functions of Ak for Ax = 450 GeV, x = 350 GeV, A = 0.83, k = 0.67 
and tan/3 = 2. 
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Figure 15: h^Q^ from turbolence as a function of A^ for A\ = 450 GeV, x = 350 GeV, 
A = 0.83, k = 0.67 and tan/5 = 2. 
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Figure 16: Plot for a (dashed line), /3/H* (dotted line), and hlfl gw from bubble collisions 
(solid line), as functions of A\ for A k = 167 GeV, x = 300 GeV, A = 0.523, k = 0.37 and 
tan/3 = 2. 
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Figure 17: tiffin from turbolence as a function of A\ for A k = 167 GeV, x = 300 GeV, 
A = 0.523, k = 0.37 and tan/3 = 2. 
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Figure 18: Plot for a (dashed line), (3/H* (dotted line), and h^Vt^ from bubble collisions, 
(solid line) as functions of A k for A x = 450 GeV, x = 350 GeV, A = 0.83, k = 0.67 and 
tan/? = 2. 
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Figure 19: tiffin from turbolence as a function of A k for A\ = 450 GeV, x = 350 GeV, 
A = 0.83, k = 0.67 and tan/? = 2. 
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Figure 20: Behaviour of the action S3/T versus T in an "extremal" region of the 
parameters-space. The horizontal dashed line indicates the value S3/T ~ 140 that must 
be reached in order that the transition takes place. 



In figs. ^ and [14] the intensity of the produced gravity waves grows rapidly as a 
function of the soft breaking mass terms. However, A\ and Ak cannot be extended beyond 
the values shown in the figures: beyond such values the condition S^TJ/T* ~ 140 is 
never satisfied and the transition does not take place. This is due to the presence of a 
potential barrier also at zero temperature, generated by trilinear terms in the tree-level 



potential: as discussed in section |I~3|, the Euclidean action S 3 (T) computed on the bubble 
solution decreases from +00 to in going from T = Td eg to T = Td es t; if the origin never 
gets destabilized - as in the cases we are considering here - one has 5*3 (T) — > const 7^ 
as T — * 0. The value of this constant is related to the heigth of the zero-temperature 
potential barrier, and it is possible that the quantity S^(T)/T never reaches the critical 
value ~ 140. This behavior is illustrated in fig. |20| []. It is clear that moving to these 
extreme regions makes /3/H* close to zero, because it is proportional to the slope of S3/T 
computed at the "crossing" S3/T ~ 140 (see fig. p0|) . By physical reasons, one cannot 
accept values of /3/H* smaller than 1: for such values the transition would be so slow 
that it doesn't complete in a Hubble time, and clearly the computations that lead to (|j) 
are not valid in such a case. Our best results give therefore values of /igf2 gw ~ 10 
peaked at a frequency, obtained from eq. |5|, of approximately 10 mHz. 



10 



t Notice that in these "extremal" cases there exists a large range of T in which S^{T)/T is low and slow- 
varying: for such cases the integral ( |l6|) must be computed numerically and leads to a little less restrictive 
condition for the temperature of the transition, however always of the form §3(T*)/T* ~ 140 — 145. 
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3.3 Gravitational waves from turbolence in the NMSSM 

We now estimate the amount of GWs produced by turbolence in the case of the NMSSM 
for the specific regions of the parameter space studied in the previous section. Substi- 
tuting in eq. (0) the values of a and /3/H* shown in figs. [12|, [T3|, [16] and [18] we obtain for 
/iQf2 gw the plots of figs. [T^, [T5|, [T7| and [0| 



We recall again that, due to the difficulties in estimating accurately the production of 
GWs from turbolence, these plots should be taken only as indicative. Still, they suggest 
that turbolence can be a really powerful source of GWs, indeed more powerful than the 
better studied mechanism of bubble collisions. In special points of the parameter space, 
we find values of h^Q^ even of order 10~ 9 , which would correspond to a signal-to-noise 
ratio of order 100 for LISA! It is clear that the production of GWs by turbolence deserves 
further studies. 

4 Conclusions 

The production of GWs at the electroweak scale is strongly dependent on the model used 
and on the parameters of the model. This is an expected consequence of the well known 
fact that the strength of the phase transition is itself strongly dependent. Our aim was 
to compare the production of GWs with the reference value given by the sensitivity of 
LISA. In the Standard Model, for the experimentally allowed values of the Higgs mass, 
there is no electroweak phase transition at all, and therefore no production of GWs. We 
have therefore turned our attention to supersymmetric extensions. For the MSSM the 
answer is again negative. The results never exceed values of order h,QQ gw ~ 10~ 16 , five 
orders of magnitudes below the sensitivity of LISA. 

More encouraging results can be obtained in the Next-to-Minimal Supersymmetric 
Standard Model, basically because in this case the transition is due to tree level cou- 
plings, rather than being radiatively induced. In this case, while in many regions of the 
parameter space the signal is still neglegible, there are also region where one can get 
hlQ gw ~ 10~ 10 from bubble nucleation and even hlQ gw ~ 10~ 9 from turbolence, with 
a peak frequency around 10 mHz. A background of this intensity would be within the 
reach of LISA. 

While regions of parameter space with such large values are, generically speaking, 
quite small, it must be observed that the condition for an intense GW production is the 
same as the condition for the generation of the baryon asymmetry at the electroweak 
scale. 
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